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e trie st’-.oil; *ti... is tne ;rcess of maintaining a set 
" refer iir- us i » * r. j of d 1st aruance 5 . Geometric 

ctsbili;;.- 1 i.. n syst . . are . .* in t. control of guided n.issi3es. 

u*c method- of geometric MZizaticn cse sir. -1? dogr^* f fre- ' n 
..tegrat y gyros to est: . . sr. inertirl refer . i dir': . Ion. 

'‘the: r. outs the gyros on a star le yl'tforni v.hich isolates the 7. 

-m mom the missile rroticr.. A second method n.o.u.ts the g/rc 5 
directly on the missile an 1 !r. effect males the missile the stable 
* latfon: . this confi mra th r angular r oti r: of the vehicle as 
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significant errors 



in maintaining reference direct! ns, hie' ir. turn loa in to error 
i r: mi . ; s i 2 e o r ie • » tn t i o n . 

Zrrors mu ile orientation ce'-eloL l e cause cf the non- 
c* u rat-, t Lvlt v . of t : ^ ? dir -u.s * : nal finite rotations and the interfering 
effects „,f Labile ang :1* r rr. ti r. o: o ro i erf r.r. ance . An error in 
^rient'-ti n can l 1 deter:, * ne d b< cor. : arii.g the angular ; ostition of 
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the missile as se jn i:i inertial axes with the angular position as 



seen in tody fixed axes. The angular velocity of the missile as 
seen in body axes can be determined by solving th* missile performance 
equations. This ~ngular velocity can be expressed in inertial ax^s 
by means of a standard transfonuati u. The angular velocity can then 
be integrated I: both cocrii ate systems ar.J the resulting angles 
com. ared. Integrating gyro gin.be 1 angles define the missile orient- 
ation for guidance. C m-. arise of gyro ~imral -r.c missile 

angular posit i n as seer: in body axes >ter: ir^s another error in 
orient^ tion. 

Tne ; e. " romance equations of the syst^.* ar' on-li: vr and 
therefor** solutions we re commuted cn a di ’ital c/ uter for particular 



sit r-ticns. The solutions showed that a systematic increase in 
error develop? d v.*ith time when the missile was sub je 't^d to ; e ri ic 
interfering torques. This effect existed reg rile us of the control 
system dynamics. Vi. :n thenmissile was subjeet-d to rar. ir.tx»r. ing 
torques the errors existed but due to the statistical nature of th* 
results no further conciusi : .-ere rr since only m single rur. • 3 
made. Further ir.vestigatlor of the system using both random and 
periodic inh -rfering torques is desirable. 
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Angle - derivatives with respect to time are denoted 
by dots over the symbol for the first two 
derivatives arid by numbers over the symbol 
for- higher derivatives. 


C d 
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CT 
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Subscript 5y bols 



Syrr.bol 


Definition 
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Body - refer to components in body nx~*o 


gim 
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Inertial - refers to component in ir.» a r!ial axes 
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Gyro unit input axis 
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All physical quantities are expressed in Gngiish ur its e <c° 
angles which are * n radians. 

xi 
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This investigate! n v;; : s conducted to determine the effects of 
angulc r base motion interference on a geometric stabilizati n system 
using body fixed gyros. Tne primary concern v.as the ability of the 
system to maintain a reference orientation in the present of errors 
produced by the effects of interfering angular metier or gyro 
performance and the effects cf thre° dimensional finite rotations 
of the base. 
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CHAPTER I 



INTRODUCTION 



Recent scientific and er gineering advances have added satellites 
and ballistic missiles to the kinds of vehicles used by man. The 
successful use of these vehicles requires a rears of cor trolling their 
motion during a part of their flight to ke?p them on a previously 
selected path. Ballistic missiles and similar vehicles ray operate 
above the earth 1 s atmosphere. The extreme altitudes require rocket 
propulsion. The operating region of these vehicles and the # •repulsion 
system used gives rise to new problems in implementing the control process 
Control implies a means of constraining and directing the vehicle to 
follow the -desired i ?th from the departure poix.t. It is the furctior of 
the control system to dire 4 the application of the r.ec es‘-a.ry constraints. 
The basic :riec*j2es sr a applicable for the er.tir* class of ballistic 
missiles * u e 1 J i n g si 1 e 11 ; t j r* . 

Rochet yropelled vehicles ex: end t:^lr * r ti re store of chemical 
fuel duri' g the first y art cf the flight. "his is folio*.; ?1 by a relat- 
ively ler ** * - rich ^f j *' o* ^ r°d fli d.t durir !/.-'» remainder f the tra : — 
ectcry. Ca x vl is applied during the y ; * ?' *- d part c p *h** flight wh=»n 
th fuel is rapidly beirg con verted into meclenl -1 'M.'upg . ~ l e c! ;’?ct 1 v*** 
is to achi'\e th** altitude ^nd velocity' necessary to ’lac** the vehicle 

* 1 c g- a 4 , th 3 menu t of **■ rus + * ?r*d* ' J 1 c n . "he 
'• 1 * r ^ ie d ? t ri r i n - ** v:i th re o v c t to its r c go. 4 t ud e 
and it*, orient iiop. ^n this res ect it is ros.rh]^ to ^rscsiat' a 



e 4 h n 



c.e on 



velccity cf x e v- 



vector with the instantaneous velocity of the vehicle. The task of the 
control system is to identify the vehicle velocity vector and maintain 
its magnitude and orientation in a prescribed manner. Most important, 
the control system must cause the vehicle to attain the predetermined 
value of the velocity vector at thrust cut-off. This means that the 
purpose of the control system is to control the magnitude and direction 
of the velocity of the vehicle at thrust cut-off. 

t 

In order to control the orientation of the vehicle's ’’viocity 
vector, the control system must (1) establish and maintain a reference 
orientation in some usable form (2) provide a means for comparing the 
actual velocity vector orientation with the reference orientation and 
(3) provide a means for altering the velocity vector to make the dev- 
iation from the reference become zero. The reference orientation need 
not be constant during the flight. It may coincide with the desired 
flight path or with any other prescribed path. It is possible to dis- 
tinguish between an inertial and an instantaneous reference orientation. 
The control system's function is to bring the velocity vector into 
coincidence with the instantaneous reference orientation. 

It is desireable from the military standpoint to eliminate any 
dependence of the vehicle on external sources of information. There 
are systems that are not self-contained, requiring some of the components 
of the control system to be located on the ground or elsewhere outside 
of the vehicle. There are other systems that require external radiation 
or external reference points to establish the reference orientation. 
Inertial guidance systems have been constructed that are completely 
independent of external sources of information.^^ These systems use 

* Superscript numerals refer to similarly numbered references in the 
Bibliography at the end of this paper. 
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measurements made vith resect to inertial space to establish the 
reference orientation and to measure deviations from the reference. 

Since the control system depends on a process of relating the velocit 
vector of the vehicle to the reference orientation (now assumed to be in 
inertial space) any uncertainty in the representation of the reference, 
or in the information that corresponds to the representation, is equivelen 
to an uncertainty in the orientation of the velocity vector. The ability 
to control the orientation of the vehicle's velocity is impaired to the 
extent that the reference orientation is uncertain. To obtain high 
accuracy in controlling the vehicle trajectory, it is therefore necessary 
to maintain the reference ori°ntation v;ith high accuracy. 

The most elementary task of the control system is to maintain the 



reference orientation v:hen this orientation is fixed in space. Geometric 



stabilisation is the ability of the system to maintain the reference 
orientation in the presence of external sources cf interference. The 
ability to track a changing reference orientation and to command the veloc; 
vector to acquire a prescribed orientation vith respect to the instant- 
aneous reference is the guidance function of the cortrol system. However 
a measure of the capability of the control system is its success in 
achieving geometrical stabilization. 

Such a stabilization system has been instrumented using three 
single degree of freedom integrating gyros mounted cn a gimbaled platform, 
with their input axes along three mutually perpendicula r directions. A 



line schematic diagram of a three axis gimbal system used for supporting 

( 3 ) 



such a platform is shown in Fig. 1-1# ' The output signal of each gyro 

is proportional to the integral of the angular velocity of the platform 
about the respective axis. The outputs of these gyros actuate 
gimbal servo drive motors. These cause the platform to rotate so as to 
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INERTIAL REFERENCE 
PACKAGE 
(irp) 

(Three single-degree- 
of freedom 
gyro units) 



[ 



INDICATED 
VERTICAL 
PACK AGE 



(ivp) 



(Two single-degree- 
of freedom 




NOTES: 1. THE BASE MOTION INTERFERENCE ISOLATION GIMBAL SYSTEM IS MADE UP OF THE 
OUTER GIMBAL SUPPORTS, THE OUTER GIMBAL, THE MIDDLE GIMBAL, THE INNER 
GIMBAL, THE ASSOCIATED DRIVE MOTORS AND THE ASSOCIATED RESOLVERS. 

2. THE ELECTRICAL POWER SUPPLIES, ELECTRONIC UNITS, COMPUTERS, CONNECTIONS, 
RACKS AND OTHER COMPONENTS NECESSARY TO COMPLETE AN INERTIAL GUIDANCE 
SYSTEM ARE NOT REPRESENTED IN THIS FIGURE. 

3. THIS ILLUSTRATION IS BASED ON FIG. 4 OF WRIGLEY, WOODBURY AND HOVCRKA (2) 
AND FIG. 6 OF DRAPER AND WOODBURY (4). 



Fig- 






-ine schematic diagram showing essential mechanical elements of inertial guidance system 
based on rotation of the inertial reference package with the indicated vertical. 
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null the gyro signals. At null indication of the gyro the platform is 
essentially restored to its original orientation in inertial space. 

Figure 1-2 is a pictorial diagram of one such instrumentation. Such a 
system may be considered as three, single axis space integrators v ' 
operating in parallel. 

The essential features of a single axis space integrator are shown 
in the block diagram of Figure 1-3. This investigation was concerned 
only with the geometric stabilization problem. For the remainder of 
this paper it will be assumed that the command signal is zero, and the 
actuating input to the system, is an interfering torcue acting on the 
controlled me nice r. The performance of a system possessing geometric 
stabilization may be described by the following steps: (l) The inter- 
fering torque acting on the controlled member causes an angular rotation 
with respect to inertial space about the gyro input axis, (2) The 
angular rotation causes a precession of the gyro with respect to the 
controlled member, producing a signal front the signal generator, which 
is proportional to the angle of precession. (3) The gyro signal is fed to 
the servo motor v:hich them applies a torque to rotate the controlled 
member back to its initial orientation, nulling the gyro out: ut signal. 
System damping is provided by compensation networks. 

The basic single axis system then, consists cf a gyro package 
mounted on a controlled member and a torque generator. The reference 
coordinate system is fixed in the controlled member. Therefore maintainirp 
the controlled member orientation maintains the reference coordinate frame 

With the reference coordinate system estaVlishei within the missile, 
the means of controlling the missile to its desired orientation within 
that frame may be considered. It is through control over the orient- 
ation of the missile, that control is ex-rcised over the vehicle's 
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supporting member 

(•«0 



OUTER SUPPORT GIMBAL 

(•••) 



BASE 

(bo) 



INTERMEDIATE GIMBAL 
DRIVE SYSTEM 



INERTIAL SPACE 
REFERENCE UNIT 

(«»rw) 



input axis 



Z-AXIS DRIVE SYSTEM 




INTERMEDIATE 
SUPPORT GIMBAL 
(••g) 



ELECTRICAL SLIP RINGS, 



».TCT! NS, HESOLVEHS ARE 'MJTrrP EH' M THIS E Id IRE. 



THE X-AXIS UNIT IS A SINGLE-AXIS INE RTIAL SPA^E R TAT V N AND ' 

the Y-axis unit is a single-axis ine htial si a t h tat: n and 

THE Z-AXIS UNIT IS A SINGLE-AXIS INERT I Al. Sf A E R ;TA IT N AND v 



MM AND SIGNAL HE 'EIVER WITH ITS INPUT AXIS ALONG X^ 
■MM AND SIGNAL RE ET.T H *!TH ITS INI UT AXIS ALONG Y. 



OMMAND SIGNAL HE 'E. IVLH • ITH ITS INPUT AXIS ALONG 



(ru)‘ 



^rw)' 



This Diagram is b^soa on Fig, 5 Of The Sherman V. Fairchild 
publication Fund Paper TTO. FF-13, Institute of the Aeronautical 
Sciences, TTew York, January 1955* (Used With Permission) 



Pictorial Diagram Of a Three-Axis Inertial Space Rotation Integrator 
3asei On Three Single-Axis Inertial Space Rotation And Comniand Sig- 
nal Receivers With Servo Drives For The Three Gimbal Axes. 



Fig. 1-2 
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Functional Pia&rak of a S/^s-lb Axis Spacb Xptsg-aat o# 



velocity vector direction. In order to function, such a control system 
must include a means of generating the desired orientation of the vehicle 
in the reference frame, a means cf measuring a ctual missile orientation, 
and a means of changing it to minimize or eliminate any deviations 
between the dcsir:d and actual orientatio: . dince tils investigation is 
concerned with stabilization, it will • e assumed that the missile has 
the ability of attaining the desired orient tio r. and the function of 
interest in the control system, is that cf mainuair ing tnat orientation 
in the : resenc^ cf inter: e ring torque. > annliei tc the fissile. 

One scheme for ; reducing control tor;u^s for missile orientation, 
uses a thrust engine to : reduce rer/torl: g torques ar out the pitch and 
yaw :>xes of the missile and exhaust ? or: cr 'he eircumfererice of the 

airframe *o central rotations : bout th : ro 12 * xis. This requires that 
the thrust enjir. lv> mounted in gimbals, giving : t 4 v:c degrees of 
rotational freak?., scout the ~ero moment ttr-st "rne * nd : system, for 
changing the engine orientation in response to command signals. Figure 1- 
is a diagram, cf ^uch a ccr figuration . 

Looking at the c err. bins tier, cf the missile control system. ?r.d the 
inertial space reference system , the ra is a means of sensing deviations 
from desired orientation, a torque generator, and e cor trolled member in 
each system. In the inertial space reference system th° gyro package 
senses deviations of the stable platform, from its desired orient* tion. 

In the missile control system, the deviation sensed, is that of 4 he 
actual missile orientation from that desired, .cth of which are specified 
relative to the stable platform. The stable platform is the controlled 
member of the inertial space referer.ee system, and the missile is the 
controlled m.emcer of the missile control system. The inertial space 
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This figure Is based on ?ig. 1-7 of "Inertial guidance" (S) a 
monograph by C.S. Draper, .7. .,'rlgley, and Sidney Lees, pub- 
lished by the instrumentation Laboratory, ’'.I.?., Cambridge, 
Massachusetts, August 1957. (Used with permission.) 



A Diagram of a ?'lsslle ,71th a Cimbaled r "hrust Motor 

Fig. 1-4 
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reference system uses servo motors as torque generators -nd the missile 
control system uses gimbalad thrust engines for that purpose, ^linirntion 
of the functional duplication of the combined systems can result in 
saving weight, and possible simplification of the ov^r-all system. 

Since inertial space reference is desired and it is the missile 
that is to be controlled, the units most reasonably eliminated are the 
stable platform with its servo drives, arid the rissi]e control system 
deviation sensing devices. This leaves the gyro package as the means 
of sensing deviations, the missile as the controlled member, and the 
missile control system as the torque gerrerstor. The principle of 
operation remains that of a space integrator, but the operating environ- 
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maintains its orientation witi. r 1 s ; e c t to inertial :;:ac^ incite of any 
base motion or the influence of lirectly act inn iut >rfe r^nce torques. 

In a s^r.se the assembly may be considered -s a bas3 motion isolation 
system. In particular it may be rgg^riod as ~n active system to isolate 
the gyro sub-assembly against angular* vibrations of the suipcrt. In a 
rocket propelled vehicle, vfners the vehicle c:ne rates intense vibr^ior 
during the propulsion interval, the gimhal support id base motion isola- 
tion system greatly reduces the effect of the vehicular vibrstior . 

This enhances the ability of the inertial reference package to maintain 
the necessary reference orientation. 

(id 



V.’ith the missile itself used ns the control*)*? d :.»rb er, the ryros 
are rigidly attached to it. In this configuration the ryros ere sub- 
jected to all the motions of the missile. 

For operation outside the earth’s atmosphere v:here aerodynamic 
forces, air turbulence and wind shear are absent, A he missile structure 
is still s ubject to tor ;ues caused 'y lech cf yrrv^ry in the inhaust 
fio;: of the thrust any* : e , r,>rents c^us-*^ i Yy th^ c' ar i i.stribution 

of th* propel];: t, ani reactions cn *he fr"? v;h*u the 4 brush 
yimra] irive v * r: cv t: ' M line of t! s . •*- * relrt'v^ to A he 

frame , Ttere ° . r° , t r * .1 v * vi) r^ 4 * on * v r -r A 1 1 a >.*hich 

are r * T i il" r 1 1' c v . ‘ s 4 ^ nr* ■. * " ■ fr-" ^ * u ■ • r ^ J . * 4 mtb j l v r : 
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Chnsihr t: r^s^cr.oO of th ^ system iascritai A i - small dislur* * r.ca , 

A s-t of in srtial -x- - rr-y ce do fir *d m- m *nc\i-”/ vi*h the ; rlr.ct *31 
be iy C/.3S of the missile before the or;s ,j t of she d i sturt--r.ee . ~ne motion 
of th** missile re lot’ ve to these ine rt 1 ? 1 ax-s r s re : 3 ter:: * r.ad . 

The equstic ns of motion for the missile 03 n be c^r/.e: i*:/ly ‘."mitten 
only ir: terms of holy axes comp onents, since it is only ir body r that 
the elements cf the inertia dyadic for the missile remain constant. If 
such a set of equations are written, they male available upon their solution, 
the components in boiy axes of the angular velocity of the missile vrith 
respect to inertial space. The corresponding components ir the inertial 
axis system will, in general, be different and can be determined ~y 
means of the standard matrix transformation for rotat'd axes systems. 

That is: 

(2-1) -'(I-BJj = V: (l-3) 3 * A (I-B) i 

Wh ere the subscripts I and B in the vector equation ( 2 - 1 ) refer to the 
components of the angular velocity vector as expressed ir. inertial axes 
and body axes respec tively. The dyadic A ^ ^ is on expression of the 
rotation matrix which is derived in Appendix A of this paper. It is made 
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up of elements involving the direction cosines relating the body 
and inertial axes set. In general these elements will depend on 
the order in which rotations about the three inertial axis occur. 
However, it is shown in Apj^endix A that for rotations which are 
small in the sen33 that the product of two angles is much less 
than one radian, he expression (2-1) becomes, in matrix form: 
( 2 - 2 ) 
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where the subscripts X^, X^, Yj , etc. i: 


! • ■ r.t i fv the 


cc 


r r a ^ ^ g of 



angular velocity (V) and the integrals of these corrpcr. ^r.ts (A) 
in inertial ar. 1 tody axes respectively. Tfr e r osier (?-'?) is 
now independent of the orle ~ cf rotation. 

The expansion of (°-2) gives for the comp on nts of angular 
velocity/ along the X-innrtial axis: 

(2-3) 



V (l-B) y = 
*1 

Integrating thi 
body has turned 



(1-3), 



- A 



(I-: 



''(I-P), 



•(i-b) v ’(:-p) 



s expression gives th' angl° thrcurh *.. - hich th-' 
about the X-ir.ertic- 1 axis: 
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V (I-B) dt 
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a (i-b) y 
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V '(I-B) it 
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The first t^rm o.. the right in (2-4) is the angle through 
which the missile has turno l about th- X body axis. Obviously if 
the fissile undergoes seme arbitrary/ rotation end th' 1 system, of 
Fig, 2-1 is successful in returning it to its initial orientation 
in body axes then the first term cn the right in (2-4) goes to zero. 
Since the gyro elements which d ' r elop the error riv'r.l are mounted 
on body axes, this is precisely vhat the syrt^- ett°;ts t > . 



This being the case, in ori r : 'or the 



,ior. -'cut the i; ^rtial 



axis to to zero, it Is require ? trv t the :ur of th** V c ~3.mdr.irg 
t-rms on the right in (2—1) *c to ’uo. It w'celd *;p^ar ‘hat this 
will nct^ in general, occur, cn } seme error in crier/ - tier, vlth reared 
to inertial s^ ace will ’eveieg Ir. the :yster. 

Th-e effect derr.cn st~rt*d above, is '?'*!/ tally g:c~ iJ ric ir 



c::' the analysis can be a:rroach°d in this light. The n f < ’- 1 has b'en 
observed in the performanc * of gyr o test ‘ .ir tables and an analvcis 

using this somewhat different approach has be^n made by Goodman and 

(5) 

Robinson, 

The equation (2-4) may be irt^r, reted from a slightly dif^er^r.t 
viewpoint which has the virtue of i.niicatir.g the nature cf the error 
in orientation which may develop. Suppose that there ar° tvo ebser- 
ers; one in the inerti-1 reference frame and one in the body reference 
frame. Each observer concentrates his attention cn his own X-axis 
and is able to record rotation of the body with resnect to inertial 
space about that axis. In addition the observer in body axes can 
cause rotation of the body with respect to inertial space about his 
X-axis. Nov: if the body is subjected to an arbitrary rotation the 
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obs^x'ver ir: inertial axes s?es rotation about his X-axis which i? 
the sum of the terms or; the right in (2-4) while the observer in 
body axes sees a rotation which is 3 ual to the first term, cn th* 
right in (2-4), Let the body now be brought to rest at some time (t^). 
Suppose now that the observer in body axes causes the body to rotate 
about his X-axis only, so as to make the ret rotation as he has obser- 
ved it zero. During this rot' tior Vrth observer, see exactly the 
same rotation. However, when th" observer in body axes has finally 
driven the body so that the net rotation as he s 2 es it is zero at 
some time (t^) , the observer in ir°rtial axes still observes a finite 
net rotation. In fact, the observer in inertial axis sees a net 
rotation: 




The advantage of the expression (2-F) becomes apparent when 
the concert of the control system of Fig. 2-1 brought into the 
discussion. The body axes set, the observer, end the means provided 
for rotating the body are exactly thi3 control system. The time 
period (0-t^) may be thought of as thn dynamic time lag in the system 
and then the rotation observed by the observer in body axis is the 
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error signal in the system* This error si r nal itself represents an 
error in orientation of the tody with respect to inertial space, but 
it is transient in nature. The system is av/are of it and eventually, 
it will be driven out. Hov/ever, there is in addition to this transient 
error, the deviation expressed in (2-5) which the system is not av/are 
of and hence, can do nothing about. The advantage of the. expression 
(2-5) is that it express ^s the error in orientation due to the effects 
of finite rotations in lependently of any transient deviations in 
orientation due to system dynamics which may exist. 

The de termination of the magnitude of the error in orientation 
discussed above, and expressed in equation (2-5) is, at least in theory, 
straightforward. The equations (2-2) taken together with the per- 
formance equations of the system, form a set of equations which specify 
the motion of the missile, v;ith respect to inertial space, in both 
the inertial ani body coordinate frames. It is only necessary to 
solve these equations for the appropriate rotations and perform, the 
subtraction indicated by (2-5), and the error in orientation is 
obtained. Unfortunately, these equations are non-linear, so th-t no 
solution in closed form for the general case is possible. However, 
for a specific system, solution by numerical methods is possible. 
Therefore, if a system which is more or less typical of these likely 
to be encountered in practice is specified, then useful information 
concerning the magnitul^ of error likely to arise in a practical 
system under the interfering effects of missile angular motion might 
be obtained. The specification of such a system and the analysis 
indicated are the purpose of the remainder of this paper. 

Integrating Gyro Errors 

Before proceding with the specification of a system for analysis 
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(2,a,8) 



one other eri'or in orientation which can develop in a systen such as 
that shown in Fig. (2-1) should be discussed. This is essentially an 
instrument error occurring in the single degree of freedom integrating 
gyro. It is of interest because its effect on the system performance 
is, in some respects, remarkably like the effect due to finite rotat- 
ions discussed above. This analysis will attem; t to distinguish between 
the two effects and effect seme comparison of their relative magnitudes. 

For purposes of discussion, the geometric effect already described 

will be referred to as the finite rotation effect^^ while this new 

(7) 

error producing effect will be referred to as kinematic drift. 

The origin of kinematic drift lies in the manner in which a 
single degree of freedom gyro performs its function. The principles 
of operation of these devices have been covered in detail elsewhere, 
so that a simole functional description of their operation will suffice. 

Fig. (2-2) shows a line schematic diagram of a single degree of 
freedom integrating gyro. As indicated in the figure, the main elements 
of the device are the soinning gyro wheel mounted in a gimbal, a vis- 
cous damper which acts between the gimbal and the case, and a signal 
generator which generates a signal proportional to the deflection of 
the gimbal with respect to the case. In operation an angular velocity 
of the case about the input axis (Fig. 2-2) causes a torque to be 
applied to the gimbal bearings which acts at right angles tc the gyro 
spin vector. This torque will cause precession of the s. inning wheel 
and its gimbal about the unit output axis . The viscous damper acting 
between the gimbal and case opposes this motion. The resulting angular 
velocity of the gimbal with respect to the case is proportional 
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line Schematic Diagram For The Single Axis Integrating Gyro Unit. 



Fig. 2-2 
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to the angular velocity of the v:Lth respect to inertial s^ace 

about the input axds. Hence in a given time interval the global vil] 



tc the displacement of the case with respect to inertial space in the 
sane interval. The signal gene r* tor produces a °igr;al pro porticnal 
to gimbal displacement and therefore its output is proportional to 



space about the input axis. 

The integrating yro v/ill perform as described above as long as 



the source of kin wV ic drift. 

The precession of the ;vro -element under the ir. fluer.ee of an 
input axis angular velocity carries the angular momentum vector cf 



momentum perpendicular to the spin reference axi of the gyro unit. 

An angular velocity of the case about the spin reference axis will 
then also cause precession of the gyroscope wheel* The output signal 
of the gyro will therefore be in error by an amount proportional to 
the product of the sine of the gimbal angle and the angular velocity 
about tne spin reference axis. 



be displaced with respect to the case by an angle which is . roportlonal 






(C-2) there 



are 





the spinning wheel with it. This produces a component of angular 



That is: 



(2-6) 
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where in (2-6) the error i3 expressed as a deviation in output 
signal# The constant is a constant of proportionality and the 
angular velocity of the case with respect to inertial space is, as 
indicated by the subscript, about the spin reference axis of the gyro, 
unit. 



In addition to the signal deviation produce! by the angular 
velocity about the spin re erence axis a further deviation is intro- 
duced due to angular acceleration of the case with respect to inertial 
space about the output axis. Angular motion of the case is transmitted 
to the gimbal and gyro assembly through the fluid of the viscous 
damper. The viscous damper applies a torcue tc the gimbal which 
resists motion of the gimbal with respect to th- 3 case. If the gyro 
wheel is assumed non-spinning so that any torque due to it does not 
act, the torque summation cn the gimbal becomes: 

( 2 - 7 ) 

I A . — A/ t \ *\ n — **C, A 

gjjn giLni (I-ca)OA d gim 

where I . is the moment of inertia of the gimbal and rvro v. r heel 
gim 

• • 

about the output axis, A(j_ ca ) is the angular acceleration of the 

case with resc>ect to inertial scace about the output axis, and c, is 
4 4 * } d 

the damping coefficient of the fluid in the viscous danger. (2-^) 
simply expresses the motion of the gim.bal assembly with respect to 
inertial space under the action of the single torque produced by the 
viscous damper. It may be rearranged tc give the first order equation 
in gimbal angular velocity: 



(2-ft) 



I . 

gm» .• 

A + A 

c d gim gim 



. .. 

' A (l-ca )CA 

C d 
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from equation (2- 0 ) it can be seen th-»t the effect )f a uniform 
angular acceleration of the case with resect to inertial space is 
to produce, in steady state, a constant angular velocity of the gimbal 
with respect to the case. (The transient period of this motion can be 
neglected in the present ca~e because the characteristic time of the 
gyro, is extremely small compared to the c*har characteristic t ? mas 
of the system.) This sli. bet*; gimbal and c~se otviou.ly a./e^r- 
ir the ch ut -‘rnsl cf the clonal genera* r a a con. tartly ircr- ~ sir.. 



si 



i t.-'Z ( her u: ifer 
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,i:n) 



> rc. a 01 . i *;.a^ : j viv.i'n 



wher* tre retie of Tmha 1 ir.-rtla to vise' us coefficient and the 



si anal g n." ra t or c * sitivit" a re lent 



~ r*. 



a single con star/: 



( 2 -?) 



D(3 •) ^ . x = k 0 / A*/ t x n , dt 
(x3-O ca 2 / (I-ca)OA 



may r.ov: v > eortbir. Z to - iv^ th- 
du : tc motion a bent Vr.° 3' in 



The e u lions (i-') cni (2- 1 ) 
total Je via tier, ch 1 r.e gyro signal 
refer rue- and th^ output axes. The resulting ^xprrssloi can be given 
the limensions cf -ng^lar displacer net by divisirg through by the 
signal generator s nusi* ivity. This is ieslreeMe sire- it is nec- 
essary to relate this deviation t. angular : i?q ?. ac^rr-^t cf the missile, 
If this is dene the total gimbal angle due to ration about all three 
axes becomes: 

(2-10) 

A gin = 3 i c u(’’ A . ) / V ‘( I-B)IA dt + D(S -lot/ 

3 > y /S/ Ue . s 

gim ,e 

v;here the second ten: on the right is the sum of (l-6) and (2-9) 
divide' by dr 'goal generator sensitivity. If the sensitivity of 
the integrating gyro unit from input ' xis angular velocity to gimbal 
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angular velocity is taken as unity ‘hen (?-20) : ..ay vrltter.: 

( 2 - 11 ) 

D(Sg). 
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S sg(Aein'.’ e ^ 
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angle in the case whs re input angular velocity *::! g:nlxl angular 

velocitv are keit -cual (l. /.. • = 1) rqr^enis the ar>*l^ 
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through ich th 1 svsten is &ttv ting tc turn the mis vile a - cut the 
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errors. Like the 


?xpr 335 : on 


(2-5 


0 th*> equation (2-11 


the effects of any 


error Ju* 


^ g , 


tvae lags in the sys‘ 



>F li’: 



eventually be drive: out. 1 he effects of system dynar ics ? re henc^ 
not indued In this expression. For a given axi: cf the rr.issii 3 , 
say the X-axi , the input axis of the gyro cor. trolling 4 d axis and 
the body axis itself are parallel. Therefore the use cf the axis 
system de pic tel in Fig. C>2) v:iI3 be iiscor.t ii.ued and hereafter the 
gyro whose input axis is parallel tc the X-axis of the missile will 
be referred to as th~ X-axis ryro. Also since the systen really acts 
to control gyro gLvbal angle and not body angle the signs in equation 
(2-12 ) will be reversed. This simply refers the error to ghr.hal 

angle instead of body angle. Equation (2-11) ray then he written for 
the X-axis in forrr sim'Iar to (2-5) as fol3ov;s: 

( 2 - 12 ) 

D(Ss)tot 



5ss(A gim’ e) 



V (I-B) v ‘ A (&im) 



= D(A - A . ) 

D gun A 
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The expressions (2-5) and (2-12) are now identical in form. 
They can be combined to give the total error in body orientation as 
it is represented by the gimbal angle. However, before expressing 
this deviation a word of caution is in order concerning the quantity 
/wgim* In succeeding discussion the reference angle for the system 
will be taken as the gyro gimbal angle, (A^ m ). As was mentioned 
above if the gyro sensitivity from input axis angular velocity to 
gimbal angular velocity is unity , then this angle is the angle 
through which the system is currently attesting tc turn the missile 



about the particular axis being controlled. In the present ca.se 
this sensitivity is unity (1C* integrating gyro) . Under these 
circumstances the total deviation of the body ax*, s from the inertial 
reference fram. after the effects of dynamic lags have been removed 
becomes the combination of (2-5) and (2-12): 

(2-13) 

D(/w - A . ) = D(A t - A_) + DCV - A . ) 

x ganrx x I B'x v E girn'x 



The expressions (2-13) along with the corresponding expressions for 
the other two axes of the system give the total current error in 
orientation in the system. In short they give expression to the 
difference in orientation between v/here the system should be going 
and where it is going. 

In the preceding paragr iphs two sources of error which can develop 

in a system such as the one shown in Fig. 2-1 have been discussed in 
some detail. The errors have been compared, and a total system error 

due to both effects has be on expressed. It has been pointed out that 
the equations of motion describing such a system arc non-linear and 
hence any attempt to determine the magnitude of these errors must 
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resort to numerical methods of solution. This rr-anr that a specific 
system must bo chosen. The choice of such a system and the solution 
of the equations of motion for the system cm the di^iVl computer 
to determine the magnitude cf the errors exrre sed i:. (2-f), (2-12), 
and (2-13) will r - token up in the remaining che; t^rs cf tnis raper. 
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The so will each b.* 1 taker up in turn, Th^ necessary r^r^o: 
equations ana the sssurr; tions required tc write them will be riven. 

In addition, Appendix P of this paper r. resents the detailed derivation 
of these performance equations. 

The Missile 

The missile is assuned to Vo a ri ’id tody. Its orientation in 
inertial srac* 1" o rrplitelv d^sorilai v y the ori of ^ s 1 1 of 

orthooor.rl ~:.°s fixv! in the ledy, Th-^se ar^ taken as pro: % Lral ax^s 
with cri-nin at t? o -issila center or r vity. relieving the -ar.ernl 
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sensitivities ir. producing roll cement in response to integrating 
gyro and rate gyro input voltages resi eetively. 

The method for producing restoring torques in r itch end yaw 
is somewhat mere com; Heated. The missile rocket rrctor is rimhaled 
tc givo it two d agrees cf rotational five Jam a! o .i the Y and I body 



ax 3 . v.,e rot - l ~r. 



.he r. tor o': out t r ? 3° ax^s ca *s^s rentier 



icr. or 



the thrust Hn^ so tlrt it o 1 : g^r -sees through th- micelle 
e ^ r t ■ r of r * v, a ^.i^. v *l^'t i on of if t r . r u a * ~1 ~ a ^ v **' ^ 3 *■ • -> 

de - s red it or*, a r * a r c ut t i . 3 r x i n t ^ 3 r j t ' ^ • 

In t: . ? s* ; 1 > ot form fe r. char, lor f r r Htin r t ■> -'‘sile 

^ , "x o 3 i ^ .*> c^d • I * y *] ► ^ * * * <a c 

3 Hf Lci- j nt iorau? in r^s'cr.o* 4 gyr* signal * c ^*v *r:r« ~h? ir.?rti* 
rea;tlor c r the ' c tor an! gndrH In rd J I*Hr.. Fm r? r.f i r rcuuir-d 
to : never, t oscillation. 4 3 u' v a pat ^ v;rv Id product r^tcr ‘ia.i a ] 
angle ir. r^s^ons-* V gvre signal* If the engine thrust l ' ^osurred 
constant, the eguatior for restoring 1 0 rogue s f_r small r.:*or glitml 
angles may he written. For : itch axis: 

(3-3) 
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Undamped natural frequency of the missile motor 
and its control system. 

- Damping ratio of the missile motor and its control 
system. 

F = Engine thrust. 
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1 * Distance from missile center of gravity to motor 
pivot point. 

pcs(ei^,A^g) » Integrating gyro signal, to notor gimtal 
angle sensitivity of the pitch control system. 

pcs(e .A ) *= Rate gyro signal, to motor gimtal angle 
r rgu J mg ° 

sensitivity of the pitch control system. 

Appendix B presents the details of the mathematics leading to Equation 3 
fcr both pitch and yaw control. 

Gyro Package 

The gy ro package contains three single degree of freedom integra'in 
gyros and three single degree of freedom rate gyres. Fig. 3-2 shows 
how these gyros are mounted with respect to body axes. Input, spin 
reference and output axes are shown for three gyros mounted on a base 
which is fixed to body axes. For the purposes of this discussion the 
orientations shown represent those of both the rate and integrating 
gyros. 

The rate gyros are assumed bo be proportional mechanisms. Angular 
velocity of the missile about the input axis produces a proportional 
voltage from the gyro at all times. The signal from those gyres is 
then of the form: 

(3-4) 

e rgu = ^rgufWje) '’x 



Whe re : 



rgu(V,e) = The input axis velocity to output voltage 
sensitivity of the rate gyro unit. 

The integrating gyros are assumed to be free of drift effects. 



Otherwise, the units are assumed to be identical to the unit described 
in Derivation Summery 2 of reference (2). Appendix B derives the 



pe rfonnr.noo enu: tion for tho Integrating gyro ur.it in ti mis of the 
gyro axes givrn in Fig. 3-3. These are then changed to give the 
equations in terms of body axes components for the orientation of 
the gyros shown in Fig. 3-3. The result for the X-body axis (roll) 
become s: 

(3-5) 
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igu 
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° £T UP 7 ~ 

sensitivity. 



Similar equations are given for the other tv:o v ody a:: 1 *: ir. A; end lx P. 

The equations (3-1 ) through (3- r ), alcn v vrith similar expressions 
to cover all body axes, co.rr letelv JefLurt the performance of the system. 
These might be combined ir a straightforward manner to give a single 
performance equation for each axis of the missile. However in 
programming for the digital computer it is much simpler to treat 
the component performance equations as given above. Therefore no 
overall performance equation for the a xes am written. 

Appendix B gives a collection of the performance equations for 
all components for all three body axes. It is convenient to make 
certain changes of the variables in these equations at this point. 

The equations of motion for the missile (3-1) about each axis 
are divided through by the moments of inertia of th^ missile about 
their respective axes. This changes the units of the equations from 
those of torque to those of angular acceleration. Trie ratios of restoring 
tcrque and interfering torque to the moment of inertia may then be 
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defined a? new variables representing the angular acceleration of 
the body due to these torques acting alone. 

The definition of angular acceleration due to restoring torque 
can nov; be carried into the equations for the torque generating systerr 
by dividing these equations through by the appro; riat^ moments of 
inertia. This puts both the body and torque generating systerr. equations 
in terms of angles and their derivatives. 

The gyro equations can b^ expressed in terms of gyro ginbal angle 
by dividing through l\ the signal generator sensitivity. This change 
also carries through to the torque generating system, ’hen all of the 



above changes have *n 
appears as a coefficient 
system. Ap;>endl>: B perf 
redefines the coefficier. 
equations along with the 
set for the problem undo 



made 4 he overall sensitivity cf the system 
in the equations of the torque generating 
orm.s the changes of variable indicated and 
ts of the equations. The resulting set of 
transformation equations (2-2) define a c ample 
r analysis. These are: 
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Where: 3 ) * Gyro gi^b^I angle to restoring angular u:c- 

' gim*‘ nr 



Mir- 



ation sensitivity of the roll control s/s tea. 

S /, ,,A ) « Angular velocity to restoring angular acceler- 

ation sensitivity of the roll control system. 

It will be noted th'-t the equations (2-2) have been mod if Lei by dropping 
the subscripts indicating body axes components. This is in agreement 
with the system equations of motion. Angles and their derivatives 
which are referred to the inertial axes set are identified by the 
subscript (i). The sensitivities defined for the roll control system 
above have identical counter parts in the pitch and yaw control systems. 

The constants required for the set of equations (3*6) must be 
specified. The complete set of thes^ constants for the system are 
given in Table III— 1. The inertia ratios were arbitrarily chosen. 

They represent what might be found in a missile of small to medium 
size such as the second stage of a two or three stage vehicle. The 
integrating gyro characteristic tines are thos* for’ th? MIT Instruments- 
tion Labor it. wy 12 gyro. 7h-s naturil frequency and damping ratio for 
the torque generating systems were arbitrarily chosen. The values 
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TABLE OF SYSTEM CONSTANTS 
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chosen are representative of what might be expected from a hydraulic 
servo in combination with the engine inertia. The system sensiti vities 
which have all been combined in the torque generating systems remained 
to be chosen. These were chosen on the basis of a linearized stability 
analysis for the system which is disc ucced below. 



In an effort to obtain some information on system* stability, the 
equations of motion for one axis were combined. All coupling terms 
from the remaining axes were assumed to be zero. The root locus of 
the characteristic equ ticn of the resulting fifth order system is 
shown in Fig. 3-3. The calculations leading to this locus are given 
in Appendix B. 

The locus of Fig. 3-3 shows two oscillatory modes. As the open 
loop sensitivity is increased from zero, the pole pair associated with 
the torque generating system moves toward the right half-plane. At 
the same time the pole pair at the origin due to the integrating gyro 
and the missile performance equation move to the left. The system 
open loop sensitivity (3 /. V 0 was chosen so as to give both 

pcSU gijn’V 

pole pairs the same closed loop damping. 

The preceding paragraphs have served to develop the performance 
equations of a simplified version of a geometrical stabilization system. 
It has been pointed out previously that the errors under investigation 
occur regardless of system dynamics. This fact can be demonstrated by 
varying the dynamics of the system just developed and solving the 
problem for the seme interfering input. Since the problem was solved 
on the digital computer the variation in dynamics can be accomplished 
by rather eli. rente ry changes in the computer program. Therefore another 
system is proposed. The body equations of the sat (3-6) v:ere retained 
but the torque generating systems and integrating gyros were assumed 



( 40 ) 




( 41 ) 



to have no dynamics. A set of equations similar to the set (3-6) is 
given for this system in Appendix B. The programming of these two 
systems for the digital computer is taken up in the next chapter. 
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CHAPTER IV 



METHOD OF SOLUTION 



Non linear differential equations are not in general amenable 

( 12 ) 

to solution in closed form. Finding the exact solution of a number 

of simultaneous, non linear differential equations is even less probable. 
However, numerical techniques have been developed for obtaining partic- 
ular solutions of differential equations, either linear or ncn linear. 

The solution so obtained will, as the name suggests , be in number form. 
For solution, the equations must be stated with numerical coefficients. 
The coefficients may be variable, but the manner in which they vary 
must be known. Certain initial conditions of the variables are also 
needed. The computation time required to solve for the dependent 
variables over a given range of the independent variable, v;ill depend 
upon the numerical technique and the accuracy required. For improved 
accuracy the increment of the independent variable between calculated 
points must, in general be decreased, thus increasing computation time. 
The numerical method used in solving the set of ncn linear 

differential equations (3-6) which are listed in the previous chapter 

( 13 ) 

is the Runge-Kutta Method with the Gill modification. ; Use of this 
method requires a known value of each variable and the first derivative 
of each, at the start of each integration step. If the differential 
equation is higher than first order the integration process must be 
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used as many times as the order of the equstion. For example if 
d 2 dx d’ 

— 9 ■ F(x, y, z, t, ~ 9 ) the value of the second derivative 

dt az 



of y is calculated at a point, say t Q> where values of x, y, their 
first derivatives, z, and t are known. The integration method will 
give the value of the first derivative of y for t * t Q - At where 
At is the integration interval* To get the value of y at the point 
t + At the first derivative, whose value was calculated for t * t 
from the previous time step w'hich started at t - At, is integrated 
using the same process. The equations which indicate the need for 
solving for lower order derivatives are called auxiliary equations. 

The eight second order differential equations mathematically describing 
the system, therefore require that the integration method be applied 
sixteen times at each time step. The transfer to inertial axes of 
rotation rates expressed in body axes give three more derivatives 
that must ve calculate:! and integrated for each time step. The equations, 
in computation form, along with their auxiliary equations, are shown 
in Appendix B. 

The basic tool used in solving the differential equations, given 

in Appendix 3 was the I. B. b . 650 magnetic drum data-procesoing 

( 14 ) 

machine.' 7 The unit used is operated by th^* Mathematics Troup in the 
M.I.T. Instrumentation Laboratory. The machine used has an internal 
storage of two thousand re isters, two tape units, floating point 
arithmetic, and index accumulators. 

The problem was programmed using a floating address assembly 

( 15 ) 

routine, called Flad, developed by the M.I.T. Instrument? tion Lab 
Math Group. The basic advantages in using this assembly routine are 
the simplifications in coding and the optmizaticn of the program to 
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reduce computation time. The differential equation routine was 
available as a subroutine requiring only entry and exit programming. 

The program could be broken down into functional parts as follows: 

1. Obtain the forcing function, for each axis. 

2. Compute the values of the derivatives and store them in 
the proper registers tc oi.t^r the l iffercntial equation subroutine. 

3. Entry and exit instruct ins needed with the subroutine. 

4. I ut the require 3 numbers in the desired form an'J read 
then out of the c.mputer. 

5. Start the cycle over again at the next time step. 

In running the ; roclem the time st^p needed to give the required 
accuracy of solution prove c to be unfortunately small. The root locus 
of the linearised system, given it* the proceeding chapter indicated 
that the postulat'd system shouli co stable. In wringing cut the 
program a sinusodal horsing functicr; o" 15 radians per second was used. 
In a linear system time increments of one twentieth of a period usually 
give accurate results. For such a time step the solution of the system 
equations diverged. Using one millisecond time steps showed that the 
solution was stable and some of the variables changed rapidly requiring 
that the time ste. be small for the digital solution tc follow the 
changes. By trial and *rvor it was found that for a time step of five 
milliseconds the commuted values of the variables agreed with the same 
computed values when one rrilliseccnd was used. At ten millisecond 
time steps the solution diverged. Therefore five millisecond time 
steps were chosen. This time step resulted in approximat.lv twenty one 
hours of computing time being required tc simulate one minute of problem 
time, using square pulse forcing functions. 
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Suggested methods of forcing the equations were by the use of 
actual acceleration data, summing sine waves, and using pulses* Use 
of data would require reading data cards at each computation point. 

This would have required the punching of at least four* hundred cards 
for each second of problem time, evaluating a sine takes 150 milli- 
seconds. For a single sin^ wave per axis almost a half a seconl would 
therefore be required each time the forcing fun^ti:r. was calculated. 

For a square pulse only additions and subtractions are required arid 
little tine is spent in the forcing function part of the program. The 
selection of square pulses usai in this invert ig? tic:: primarily 

det^minid by the computation tire required. 

The system was programmed to determine the response to three 
forcing functions. 

1. An initial d in lac ur.enl of the rissiie about all three 

axes. 

2. Feriodic rectangular pulses of angular acceleration 
applied to <0 1 three ax 3. 

3. Random rectangular pulses of angjlar acceleration acting 
on all three axes. 

The initial displacement of the missile about all three axes was 
usid to determine the transient behavior of the system. A similar 
initial angular displac * .ent condition could be realized in an actual 
missile. For a multistage missile the sej aration process would be 
expected to cause some rotation of the next stage, b’hen the propulsion 
system fired the control system would see the displacement as an initial 
error in orientation , 

Rectangular pulses were used for periodic forcing of the system 
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primarily da? to the speed with which th -y could bo go no rated, as 
noted earlier. Use of these pulses added no appreciable commutation 
time over that required for the initial displacement, case. These 
pulses had a magnitude of four tenths radian per seconl squared. 

For each axis the pulse duration was a tenth of a seconl with two 
tenths of a seconl betws'n the -n ! of cr.- juls*- anl the start of 
an^ f h*r. The pulses were phased so only one axis at a t*m° was being 

forced, idler thi pulse was taker of. J of one axis and applied to 

another, its sign was changed. This gave the juried of the forc'ng 
function's fundarentyl frequency as six ter ths of a second. 

The random pulses of angular accelereti n -ere limited in 
magnituie from zero to ninety mil.*! ir^dians per second squared, in 
increments of ten. Pulse duration was limited to five hundred milli- 
seconds in ir.crenei ts of five milliseconds. The sign, magnitude and 
duration of each pulse was detemired using digits fzcm registers 
containing calculated values of seme variable. These registers 
contained values in floating point form. ho digits less than two 

digits to the right of th? first significant figure w?re used. These 

digits varied at every tim ? step. The digits used to determine the 
forcing function for one axis cam** from variables primarily influenced 
by motions about the other axes. A single digit determined the pulse 
magnitude. To determine the sign of the pulse, five v/as subtracted 
from a single digit and the sign of the difference was used. The length 
of the pulse ir. five millisecond increments was date mined using tvro 
digits. These procedures give a forcing function on each axis which 
is independent of the mode of forcing on the other two axes. The 
coupling between axes given by Suler's equations and the gyro equations 
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does not however make the angular motion -bout any one axis independent 
of the motions about the other two axes. 

Programming for the further simplified systems which were mentioned 
in Chapter III was accomplished by removing parts of the original 
program. 
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CHAPTER V 



discussion cf results 

This chapter discusses the results obtained from the solution of 

the system equations on the digital computer. The discussion is concerned 

primarily with the angular differences developed in Chapter II, These 

differences represent errors in indication of body orientation. The 

difference l e tween inert ial angle and body angle is the error in 

orientation due tc finite rotations. The difference between body angle 

and gyro gimbal angle is the error in orientation caused by gy ro 

characteristics . The sum of these differences ( A T -A . ) is the net 

I gun 

error in orientation due tc these effects. These differences do not 
include the deviation in orientation re suiting from dyna . ic lags in 
the system. 

Transient Re sponse 

It was de s ire able that some examination cf systen stability be made. 
In order to do this the response of the system to an initial angular 
displacement was computed, ho interfering angular accelerations were 
applied during this response. 

Figs. 5-1, 5-2, and 5-3 show the response of the system in restoring 
from an initial angular displacement. In these figures the gimbal angle 
is taken as the reference variable. This was done because the gimbal 
angle is the variable which the system attempts to null. The figures 
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for the pitch and yaw axes show a damped oscillatory respcr.se. The 
response in roll was only run until the gimbal angler, on the other 
two axes were essentially zero, Stability of the roll axis system 
was not a problem since proportional control was used in th ? torque 
generator and the loop sensitivity was relatively low. 

The root locus of the linearized pitch control system was discussed 
in Chapter III and shown in Fig. 3-3. The close loop gain gave a 
damping ratio of 0.445 for the two oscillatory modes cf the system. 

The control system* in pitch and yaw is a fifth order system in angular 
velocity as shown in Chanter III. The system is also non linear, but 
some feeling for system characteristics may be obtained by comparing 
it with a linear second order system. Charter 19 of Ref. 6 contains 
curves of damping r tic versus transient peak ratio for linear second 
order systems, entering there curves with the transient peak ratio 
obtained from the system response gives a damping ratio. A damped 
perio. of oscillation may also be obtained from the transient response. 

An undamped natural period may be calculated using the damping ratio. 

The transient responses shown in Figs. 5-1 and 5-2 are not those 
of a second ,rdar system. Therefore values o^ damping ratio and 
undamped natural frequency based on second orler res onse characteristics 
are meaningful only as a rough measure of the system characteristics. 

From Fig. 5-1, the yaw axis response indicates a damping ratio cf 0.48 
and an undamped natural frequency of 1.19 cycles per second. From the 
pitch axis response, Fig. 5>-2, a damping ratio cf 0.46 and an undamped 
natural frequency of 0.95 cycles per second were obtained. In Fig. 5-3 
the slight recurvatuxe of the gimbal angle curve indicates that the 
roll control system is rather heavily damped. The transient response 



( 53 ) 



indicates that the system is dynamically stable on all axes. 

The differences between inertial, body, and gimbal angles are 
plotted in Figs. 5-1, 5-2, and 5-3. No interfering angular accelerations 
were applied to the system during the transient response. The existence 
of these differences during the transient response indicates that any 
angular motion of the missile may cause errors in orientation. The 
absence of any means in the system of sensing these differences is 
shown by their constant value once the gyro giirbal angles have reached 
their steady state null position. 

Information on the source of the differences can be obtained by 
considering th^ orientation of the gyro axes relative to the missile. 

To aid in comparison, the body axes orientation of the gyro s are 
tabulated below. 



| Gyro Axes 



Rotation Sensed 
Roll 
Pitch 
Yaw 



IA 

X 

Y 

Z 



SRA 



CA 



Y 
X 

Y 



Z 

-Z 

-X 



J 



Table 5-1 

INTEGRATING GYRO ORI STATION 

In the equations (2-10) for the integrating gyro, the terms 
introducing differences between body and gimbai angle are angular 
accelerations of the missile about the gyro output axis and the product 
of gyro gimbal angle and angular velocity of the missile about the 
gyro spin reference axis. 

Comparing the body to gimbal angle difference curves for the roll 
and yaw axes shows a marked similarity in the oscillatory character of 
the differences. Both gyros are oriented with their spin reference 
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axes along the pitch (Y) axis of the missile. This would indicate 
that the dominant cause of this difference was the result of angular 
velocity about the spin reference axis. Th° greater magnitude of the 
error in the roll axis plot can be accounted for by the fact that the 
ginbal angle remains larger for the roll axis than for the yaw axis. 

The relatively small effect or; the body to eimbal ancle difference 
resulting from angular accelerations about the output axis is indicated 
in the pitch axis plot. The pitch axis gyro has its spin reference 
axis along the roll axis. The slow rate of change of gimbnl angle for 
the roll axis gyro shews that the roll rate of the missile is low. 

This indicates that the titer: gyro is little affected by angular 
rotation about its si in re fe re nee axis. The output axis of the pitch 
gyro is along the yaw (Z) axis. The curvature cf the Z-axis gyro 
gim.bal angle curve indicates the magnitude of the angular acceleration 
of the missile about tnat axis. It will be noted that the peaks of 
the oscillations in tne pitch axis body to ginbal angle difference 
curve occurs when the Z-axis ginbal angle curve changes direction. 

In this transient response the angular acceleration of the body about 
a gyro output axis appears to result Ir. little error. However, no 
general statement sa vms warranted for the case of arbitrary missile 
motion. 

The inertial to body an gle difference curves are almost identical 
for the pitch and yaw axes, while this difference is almost nonexistent 
in the roll axis plot. This difference is the result of finite rota- 
tions as was shown in Chapter II. For the roll axis the difference 
equation is: 



( 5 - 1 ) 
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For t!v other axes the difference is likewise the difference between 



the products of „ ody angle an: angular velocity components associated 
v. r ith the other two awes. 

Comparison of Figs. 5-1 and f>-2 shows that the gimba] angle curves 
are almost identical. This moans that at any time during the response 
the body angler and rotation rates are nearly identical for the Y and 
Z-axes. Tt e two tones ir. the difference 3 qua tier above therefore suit. 
to almost zero. "his also makes the pitch ar.d yaw axes differences 
almost identical witn time as shown in the figures. 

The sy s ter. res; onse to an initial angular displacement could hate 
a parallel in an actual missile. A similar event occur in the 

case cf a second or later stage of a multiple stag-* missile. In the 
interval between separation and igr.iticn of the next stage, a transient 
motion of the missile may be expected. This transient, motion imposes 
the initial condition on the central system wh-r. the thrust engine 
fires and control is coimmenced. If the angular velocities are small 
and car. le neglected, this woul i establish initial conditions similar 
to those used for the transient res; onse. The body vculd have a 
deviation in oriental ‘cn with zero rotation rates. For the response 
calculated the difference between inertial angle =nd body angle was 
approximately cne per cent of the initial displacement, for all three 
axes. This difference is the total orientation error developed during 
the response . 

The information obtained from this calculation of the transient 



response of the system is as follows: 

1. The control system is dynamically stable, 

2. Rotations of the missile produce errors in orientation 
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due to the effects of three dimer. s^onal finite rot' ti^ns 1 the 
characteristics of the integrating gyro. 

3* These errors in orientation are net sensed by the 
system and exist in steady state, 

4, For the response computed the total steady state error 
in orientation was approximately one per cent of the initial displace- 
ment . 

Periodic Forc ing 

The system response to a eriodic forcing function in interfering 
angular acceleration was obtained. For reasons already discussed in 
Chapter IV all of the forcing functions used in tnis analysis w^re 
composed of combinations of square pulses. For the periodic forcing 

function trie se :ue:. *e of :ulses shown in Fins. 5-4, 5-6, and 5-8 was 

arbitrarily selectee. On a given axis of the system these pulse: 
are alternately positive and negative with a period of si:-: hundred 
milliseconds, nach pulse nas a duration cf one hundred milliseconds 
and a magnitude of four hunired mill i radiar s per second squared. The 
phasing of the pulses cn trie various ax-s was arbitrarily chosen so 
th~t the yaw axis w^s tv:c tuirds of a period behind the pitch axis 

and the roll axis was two thirds of a period behind the yaw axis. 

The periodic forcing function just described is equivalent to a 
linear vibrational acceleration of the missile of four feet per second 
squared at a distance of ten feet from the center cf gravity. In the 
absence of system restoring torque this acceleration acting for the 
duration of the pulse would produce a displacement of about a quarter 
of an inch at ten feet. 

The response of the system to the periodic forcing function just 
described is plotted in Figs. 5-4 through 5-9 for all axes. These 
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Yaw Ax»s Uiiyou 



pjLots shov; the forcing fund J.c-r , the ro gimbal angl? ? rd the thi*ee 
angular differences for each axis for the first ten seconds of system 
operation. At the end of ter*, seconds the system had reached steady 
state on all axes so that further computation was unnecessary. 

The plot of the gyre girrtal angle is essentially the same for 
all axes. After an initial transient period the gimbal am-le settled 
down to a steady oscillation about its r ull position which appears to 
be almost sinusoidal. However, comparison of these t lots with those 
of Figs. 5—1 through P-3 shows that the gimbal angle response is 
actually made up of a continuous serifs of exponential curves similar 
to those obtained for the transient res pen e. 

The behavior of tne gimsal angle during the initial transient 
period is consistent with the results of the transient analysis 
discussed previously. 2ach axis develops some dynamic error during 
the transient period which is eventually damped out. In the cas ■ of 
the pitch and yaw axes this daiqing requires approximately two seconds 
while for the roll axis approximately six seconds are re paired. 

The effects, of coupling among axes in the system can te seen dr 
the plots of gimbal angle for the piten and roll axes. Sach of these 
axes is initially unforced and pet in ^ach case some small gimbal 
angle develops immediately due to the motion of the missile about the 
yaw axis. 

The angular difference of Fig. 5-4 through 5-9 all develop in 
a fairly consistent pattern in spite of the variation in dynamics 
between axes of the system. Some general discussion of the basic 
patterns is warranted before considering the results axis by axis. 

The difference betv/een body angle and gyro gimbal angle is, in 
each case, an oscillatory curve which slowly drifts a way from zero. 
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In steady state this is a combii ation of a steady oscillation about 
zero and a constant drift rate. Tnis can be explained in terms of 
the two error producing effects discussed in connection with the 
integrating gyro. 

In steady state the components of angular acceleration of the 
missile are symmetric about zero in the present case. Therefore for 
a full cycle the net angular accelera tioi of the missile about a 
given axis is zero. The difference between body and gimbal angle may 
oscillate due to the effect of angular acceleration about a gyro 
output axis but for a full cycle the net difference produced by its 
effect is zero, Any drift ir the ’ lot of body to gimbal angle 
difference must them be due to the effect of angular velocity of the 
missile about the gyro spin reference axis. 

The product of gyro gimbal angle and missile angular velocity 
about the gyro spin reference axis v/ill not in general ir.tegiate to 
zero over a full cycle. Therefore over a full cycle a net difference 
between tody and gyro gimbal angles will usually exist. In steady 
state the different due to this effect will be the same for each 
cycle. This produces the linear increase upon which the oscillatory 
effects are su x > c .e s . 

The carves for the difference between inertial and body angles 
are characterized by a small oscillation, the mean of which is linear 
with a finite slope in steady state. This linear divergence is to be 
expected for the case of periodic forcing. A net dif f > re nee developed 
during a given cycle of the forcing function would also be developed 
in any other cycle in steady state. Therefore there would be a linear 
accumulation of these differences with time such as is shown in the 
curves. 
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The sum of the differences just discussed is the difference 
between th3 inertial ar.^l? and the gyro gimbal angle. This difference 
is also plotted for each axis. It represents the total error between 
the reference orientation and the orientation to which the system is 
currently attempting to drive th^ missile iue to all of the effects 
under investigation. 

The angilar differences for each axis at t .e end of ten seconds 
are given in Table 5— 2. Since all differences steadily increase in 
magnitude these are the maximum errors in crie- t* tion which the system 
experiences for the ten second run. For those niff re:c?s which a re 
oscill' tory the numbers given are tr.e mean error at ter. seconds. 

Table 5-2 also giv-^s th d ratio of the angular diff-? ranees to ‘.he 
maximum amplitude of tne missile angular oscillation ir. steady st te. 
3ince the gyro sensitivity from bo iy angle fc v g/rc girr.bal '-ngle is 
unity, th* . amplitude of gyro gimbal an 'ie oscillation can be 

used in forming ratios. However, it is more .meaningful : o relate 
errors in orientation to body angular motion. It shouii be n^ted that 
the magnitude of the body angle is i.. t that of the gimbal angle due 
to the errors. However, the amplitude of the body angle oscillations 
are the same as the amplitude of the gi;, r 1 angle oscillat i cnr> . 

The ten seconds of running time for which the results of periodic 
forcing have been ietermined is short com, a red to the firing time of 
most missiles. A good approximat ion to missile firing time would l ^ 
ore hundred seconds which would require a computer rur. ten times as 
long as the one that was mad?. However, once the system has reached 
steady state operation and remained there long enough to establish 
the trend of the errors further commutation is unnecessary. The error 
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curves, being linear, luay be extra; >1' t id to any desir'd value. In 
the present case this extrapolation has be^n carried out an 1 the 
resulting errors ; n orient tion are sho\/n in Table F— 2 a lor. g with 
the values for ten seconds. The values of t.he differences at one 
hundred seconds are given in both radians and degrees for convenience. 

The following infonaati- r. was o' r tain ? 1 from the response of the 
system to a periodic interfering angular a :celeraticn: 

1. Errors ir. orientation develo.ai due to the effects cf 
thr-e 3 dimensional finite rotations and th? chs racteristlcs of the 
integrating gyro. 



2. These errors v/ 3 ro systematic, and ir. steady state, their 
mean increased linear]., with time. 

3. The total error in orientation developed at the end of 
10 seconds was 0. A 1 r.i 1 lir: dia.as in roll, 1.0? mil] irr iians in pitch, 
and 1.49 miiliradiaus ir ya. . These are respectively 341, 2^, and 

401 of the amplitude cf l issile angular motion about the axis in cuestic 

4. Jince the increa.se in error vr, a linear, extrapolation tc 
100 seconds was ossibie. This increases the val^-s above bv a factor 



of ten. 



Random Forcing 

In order to simulate ir.te rfering angular accelerations as 
realistically as possible within the limitations already imposed by 
computing time, the response of the system tc a series of random 
pulses was determined . Since no steady state could be expected for 
such a forcing function the problem was run for forty seconds of 
real time. This required approximately ten hours of computing time. 

The random forcing function was r.ade random with respect to 
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three of tha pulse charact*; rlsbics . The sign f t..- ,/als?, Its 
magnitude, and its duration Wire independently chosen in a random 
fashion. This was done automatically in the computer. The sign of 
the pulse was chosen first. The mechanism ~y which this was done v;as 
such that positive and negative signs each had a ; robability of 
occurrence of one half in a riven selection. The magnitude of the 
pulses was limited to the values from 0 to 90 illiradians per second 
squared, in steps of ten millirad'uns p~- second squared, nach 
value had a probability of occurrence of cne tenth in a given select- 
ion. The duration of the pulses was select ;1 at random in the range 
from 5 to 500 milliseconds with selection be i n a limits I to multiples 
of five milliseconds. The forcing functicn obtained in this manner 
is shown for eac.n amis in rigs. 5-10 through 5-12. It should be noted 
that trie re is n c correlation between the forcing functions for the 
various axes since selection or. each axis was completiy i:;i\endent 
of the other two. 

The response of the system to the random forcing function just 
described is snown in Figs. 5^10 through r-13. 

The system response as shown by gyro gimbal angle is about what 
could be expected in view of the previous results. The gimbal angle 
oscillates randomly in essentially the same manner as the forcing 
function. There is no steady state condition apparent in the response. 
None was expected since the forcing function is random. In general 
the behavior of the gimbal angle is consistent with the previous 
observations on system dynamic stability. 

The angular differences have the same general character as those 
observed in the periodic case. The body to gimbal angle difference 
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Total Error in Oriintation with FaROon Forcing Function 
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curve is oscillatory, while the inertial to body angle difference 
is relatively smooth. However, there is no systematic divergence 
of these angular differences for this run . In this resect the 
response was different from the results for periodic forcing, 

3ir.ce the response to random forcing is essentially statistical 
in nature generalized statements can not be made on the basis of a 
single run. The angular difference curves shown are just one set of 
an infinite number that couud occur as a result of this type of 
forcing. The fact that the angular differences obtained during one 
run do not diverge, does not imply that th?y would net diverge in 
another run . However, the sir.rle run shown does demonstrate that 
errors in orientation do exist ever, with random forcing. The 
inte rf ering angular acceleration inputs nee; r ot be periodic for 



these orientational errors to develop; and these are the only conc- 
lusions th t can reasonably be drawn from the re salt" cf a single 
run. ho re runs would permit a statistical approach to- the problem. 
More generalized results could then be obtained. 

A run was made to show that orientational errors will exist 
regardless of control system* dynamics. The system dynamics were 
altered by assuming proportional response in the integrating gyros 
and the torque generating systems. The reduced equations are shov .71 
in Appendix 3. 

This simplified system was forced by the same periodic forcing 
function used previously. The response to this forcing is shown in 
Figs. 5-14 through 5-16. 

It will be noted that the plots of this response show missile 
body angle as the reference variable. For this simplified system 
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Response in Va\n fop System with Reduced Dynamics 



the gyro angle is identical to the body angle. Therefore the only 
difference that can occur is that between inertial and body angles 
The angular difference curves from inertial to body angle indicate 



the error in orientation develops in the same manner as for the 
complete system. The error is slightly less in this case since 
the control system response is faster. Therefore the missile will 
experience smaller angular deviations and lower angular velocities 
for the same forcing function. This accounts for the reduction in 
the net error in orientation due to finite rotation effects. In 
any case the error does exist regardless cf control system dynamics. 

The results obtained in this investigation prove the existence 
of errors in orientation due to thr a e dimensional finite rotations 
and integrating gyro characteristics. Some inform ation was obtained 
for the cases of per iodic forcing, random forcing, and transient 
response from an initial angular displacement. 

For the case of periodic forcing, further investigation is 
needed to determine the effects of the frequency and magnitude of the 
forcing function. In the case of random forcing more runs will be 
required to permit a statical evaluation of the errors developed. 

If very many runs are to be made, a faster me h tod of computation is 
highly desirable. This can be accomplished through the use of a 
faster computer, or perhaps by changing the method of computation. 
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CHAPTER VI 



ai;c:u.5ic::s a:.d r<:c. 



)A7I0KS 



The preceding analysis has shown that 
may develop in th d : n ^ 4 rical siatili ’' tic 
gyros. These errors are iistinet from th*^ 
which such a system would experience due to 
was demonstrated that urvdsr certain cor.diti 
ation systematically increase v.'ith time. A 
errors in orient tior. develop regardless of 
Errors ir. orient tier developed in the 
angular vibrations. In restoring frc-ir. an i 
the system developed an error of about .ne 
displacement.. This error remained when the 
Those errors dev_*lo 4 .eu due tc th- 3 coupling 
between axes of she system. 



errors ir orientation 
n system using body fixed 
\ ‘-J v : a 4 i. r.s in orientation 
system dynamic lags. It 
ons the errors in orient- 
ls. , it was shewn that 
system dynamics . 
absence cf interfering 
nitial angular displacement 
: e r cent of the initial 
system had come to rest, 
of missile angilar motion 



then the system was forced with a periodic interfering angular 
acceleration or. all res, t re re was a systematic increase in orientation 
errors v.ith time. After ten seconds trie total error in orientation in 
yaw was forty per c^nt of the r. aximum amplitude of the angular oscillation 
of the missile in steady state. Corresponding figures fer the pitch 
and roll axes v:ere twenty seven er cent and thirty four per cent 
respectively. These errors increase linearly with time. Therefore at 
one hundred seconds the error in yaw becomes four hundred per cent of 



( r 8) 



the n*aximum amplitude of the angular oscillation of the missile. 

The same kind of errors in orientation developed when the 
interfering angular accelerations w«re random, however, due tc the 
statistical nature of the forcing function no conclusions could be 
made on the basis of a single r in concerning the r.ar.mr in which the 
errors will accumul' t“ a i th time. 



rurth-fr ‘ nv *stig -.ilon cf til** : r hlix is still required. The 
investigation of ,eric*ic fencing should be extend'd tc dote ranine 
the effects of the frequency and magnitude cf the interfering anguia 
needier tions. Zn^ugh results sh:-ul 1 * e obtained u-ln* a random 



forcing f 11 ctlon to : 


emit a statistical 


r- r. 


n- lysis of the errors. 


These additional inv- 


• =: t i ’ u t i on s would t_e 


O' 


eatly facilitated if a 


sign! f leant reduet lor 


: in con put in ; t ? r* 


CO J 


ild achieved . 
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ArFSMDIX A 



TH.3 TRAKSF.SR CF VVXTcR C^I I C 11.73 HITV73 5N XT AT *2 D AX33 5TT3 



Fig. A-l show *3 the inertial and body axes sets described in Chapter 
II and an arbitrary vector V which is fixed in the body axes set. The 
roblen is to \xj the coi.yor.ents of the t< r in the inertial 

a,.es set ulvm the body axes coirron-nts * known and the body axes, 
have assumed an arbitrary orientation. in general, an arbitrary 
order, tat lor- of th^ body axes relative tr. the inertial a ■ ?s can be 
duplicated by, t u-st, three succe. sive fin’d ? rotations. Ir. th- 
analysis belcv. three successive rotations of the body, one about each 
of the inertial axes, are assumed. After e = ch natation the conjcnents 
cf the vector V in inertial a.e s are derived. After three rotations 
the results for th.e Individ r* 1 rotations are combined to give the 
components of the vector in inertial axes in ter»s cf th° three 
angles of rotation and the comLorur.ts in body ax^s. The resulting 
transformation ecu- tion is not unique. There is a different trans- 
formation for eacr. possible order of rotations. However, a simplification 
is made at the end which gets around this difficulty. In the derivation 
telov: the order of rotations is about the X,T, and Z-axes in that 
order. 



Let the body axes of Fig. A-l be 



axis through the angle A y as shown i 

h 



rotated about the X inertial 
Fig. A-2. From the figure it 



(SO) 



can be soon that the compor. nts of the vector V in inertial axes 



can now bo written: 

(A-l) 

V = V 
*1 *B 



V » V cos A - V sir. A 



y I y B "I 



Z B y T 



I * V sin A + V cos A 
Z I y B X I Z B X I 



There is consider' hie simplification afforded by wril 
equations (A-l) in matrix form. 

(A-2) 

0 0 



V ' 1 1 



X I 1 
V = ! 0 

y I ' I 



..rt — >/•. 

X T X-r 

1 J. 



V 



V 



i V 



0 BA 






I 



'B 



’vhere the notation 5A , and CA r?fs:s to the sir.e 

\ y T 

the angle of rotat i n about tie i-axis. 

h ov; let a new set of body axes (B } he defined : 
are coincident wits the inertial axes of Fig. A-2. I 
transformation (A-2) always applies the sets 

(B) and (B*)» That is: 

(A-3) 

I 1 0 0 V 



Vx B 



*B 



I | j 

V' : = 0 CA -BA V ! 

y 3 *1 *1 y 3 



1 V* 

i z r 



3A CA V ! 
*1 *1 Z B 



l B J L 

If the body is now rotated about the Y inertial 



ir.g the set of 



and cosine of 

such that they 
.cte that the 
of : oiy axes 



axis through 
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A then from Fig, A -3 it can be seen that: 



(A-4) 





l"v “ 
X I 




CA 

y I 


0 


3 A 

y I 


*b 




V 

y l 




0 


1 


0 


V' 1 
y B | 




V 

Z I 

_ 




-3A 

y I 

i. — 


0 


CA 

y T 


V» 

3 J 


and def 


ning the 


body ax 


es set 


(3") in 


exactly 


(3') gives: 












(A -5) 














: 


R..n 

' ^ | 


i CA 
: y I 


C 


3 A 

y I 


V 

*B 




V ' ! 

y n 1 

LJ , 


1 


0 

1 


1 


0 


^ t 


: 


V ' , 

L 


! ; 


I y i 


0 


CA 

y I 


y t 

JBJ 


Rc 


>ta ting 


nov: about 


Z T through A 
I z T 


gives: 


(A-6) 


r v * 

i Xl 












i 


| 


' CA 

Z I 


-5A 

Z I 


0 


v* * i 

X B j 




V 

y I 


i 


3A 

2 I 


CA 

Z I 


0 


i 

V' 1 } 

7b 1 


j 


V 

L \ 


1 


0 


0 


1 


V f » i 

Z BJ 



Replacing the matrices above by symbols of the form (V^) allows 
equation (A-6) to be v.ritten: 

(A-6) 



(Vj) = (A z ) (V”) 

where the meaning of the notation becomes evident by comparing the 
two forms of the equation (A-6). 

If (A-3) is substituted into (A-5) and this equation in turn is 
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substituted into (A-6) the result is: 

(A-7) 

(Vj) - (A 2 ) <A y ) (A x ) (V B ) 

If the matrix multiplication (A ) (A ) (A ) is carried out the resulting 

x y z 

equation for the inertia] components in terms of body components 
be c ome s : 



(A-*) 



1 V 




CA CA 


oA C).-\ oA — S-rV 


CA SA CA +3A 3A 1 h 


I *1 

1 


. 


2 i y i 

I 


Z I *’ I *1 Z I 


Z I y I X I Z I *1 x 3 


V 


1 J 


1 3A CA 


SA 3 A 3 A +CA CA 


^>A on CA -CA o A V 


71 

1 


! i 


! *1 ^ 


y i *1 z i z i *i 


'“I y I X I "i X I y B 


1 V 


[ i 

i 


t 

*-oA 


C A On 


CA CA V 


. Z I 


! 


L iT 


y I -'-j 


v X 7 

"I X I , Z 3; 

L- Jt 


If 


the rotate cn 


3 a boat the inertial ax 


? 3 ar Kept smn I j , so that 



the product of th * sin -> of two of the angles can be neglected this can 
be written: 

(A-9) 



l V ' 


r 1 -A A ~ ! 


f v 1 


*1 


“I y I | 


x b| 


V 


= A 1 -A 


V , 


y I 


1 I X I 




y B 


V 


-A A 1 




V 


“I 


r 

i 


i 

i 


Z R 

L 1 



It will be noted that the restriction imposed on the size of the 
angles in (A-?) is more severe than that imposed by the us cal small 
angle assumption. Therefore replacing the sine of the angle by the 
angle and the cosine by unity is justified. 

It was noted above that the transformation given by (A-^) is not 
unique. In fact there are six different matrices similar to the one 
in (A-3) corresponding to the six orders of rotation possible about 
the inertial axes. However, under the small angle assumption made 



these all reduce to (A-">) . The equation (A-9) and (?.-?) ere identical. 
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Initial OhjBhtation of Bcof Ax&j JLn£RTial AxFS^ aa& fecTo/? V 

f't A ' i 
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0*!EnTat/om A(?T£R Rotat/oa/ Asoor X £v ER tiAl Ax'S 

Fif. A-£ 
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Orientation After Rotation Aqoot Y Inertial A\is 

F'f A “3 
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Final Orientation after. Rotation About zZSn£/?T/c Axes 

F '% A ‘? 
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AFFSkDIX B 



SYoTi PERFCH:-:AKC3 s.u/.tichs 



The performance equations for the system of Fig, 2-1 are developed 
belo'A in the order of their discussion in Chapter III. 

The missile is taken as a rigid body. Since translational motion 
is of no interest in this analysis , the center cf mass of the missile 
is assured fixed in inertial space. The • issile equations of motion 
are then obtained by summing moments about the missile center of 
gravity, expressing components in body axes gives Zuler's equations: 



(3-1) 



v; + (i -i ) v; w = !•: 

x z y z y > 

w + (i -i ) v: v; = h 

y x z x z j 

V. + (I -I ) W V = M 

z y x' y x 2 



where in the present case the sum of the a 
above includes the restoring torque ap; lie 
control system and the interfering torques 
vibra tier . Thus : 



polled moments on the right 
d to the missile by the 
due to missile angular 



(B-2) 



M 



x 




{•; 

y 





+ 



‘"‘(int) 

x 



+ i, 

+ K 



(int) 



y 



(int h 

(3 9) 



The remainder of the system performance equations are concerned 
with the determination of the restoring torque (M^) as a function of 
missile motion. 

In the missile configuration being considered restoring torques 
for the pitch and yaw axes (y and z axes) are produced by rotating 
the thrust line of the missile engine away from the missile center 
of gravity. Restoring torque is roll is produced by two nozzles on 
the circumference of the missile at the opposite ends of a diameter. 
The components of restoring torque about the pitch and yaw axes may 



be written: 

(B-3) 



K = T 1 A 
m m 

y y 



m 



T 1 A 



in 



z z 

Where: T « the thrust of the missile motor, assumed constant. 

1 = distance from the missile center of gravity to the 
motor pivot point. 

= the angle through which the rotor is rotated about 
the pitch or yaw axis (Assumed small) 

Similarly the r? storing torque in roll becomes: 

(B-4) 

M - T, \ d 

m x ( re ) 



V.Tiere: 7, x = the thrust of the roll engine. 

(re) 

d = the diameter of the missile. 

The system which positions the missile motor in pitch and yaw 
is assumed to be a second order system which receives gyro signals 
as input and positions the missile motor about the appropriate body 
axis. 
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The performance equations may be written: 



(B-5) 



+ 2(DK) 



pcs 



pcs 



pcs 



3 



III 



v; 



, 2(DR) 

+ v v 



yes 



ves 



yes 



+ 1 



S pcs(e i ^ u ;A) e (i«u) y 

^Dcs(e ;A) e (rgu) 
y rgu ° y 



"ye s(e igj ;A) e (i^u). 



+ S 



.vrere : 



ycs(e ; A ) (rgu) 

" rgu J ° z 



(Dl)^ = 7 itch ccr.trol system tamping ratio. 



=* } itch control system natural frequency 



ICS 



3 / , \ = r itch control svstex ser.sitivitv from 

ccs(e. ;A) 

lgu’ 

integrating gyre signal to motor angle. 

^nesfe :A) = 7 itch c cr.trcl s stem ssnsitivitv from 
rra 3 

rate gyro signal to motor angle. 

Similar .tafinlti jr o hold for the yaw control system. 

loll engine thrust is assumed to be directly pro* ort ior,al to 
gyro signals. Tnat is: 

(3-6) 

v e . + S 



*(re) rcs(e. ;?) igu 
i?u x 



rcs(e ;T)'rgu 

rgu x 

Zquatio: s (3-3) and S-f may a combine 1 to give the restoring torque; 
ir. .r tee. and yaw as a function of gyro signal. 2quati<>r.s(?-4) and 
(E-6) may be similarly combined. Then: 

(3-7) 

= °(rcs)(e. ;X) e igu + °(rcs)(e ;M) e rgu 
x igu w x v rgu * x 
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2 



V 

n 

— P cs 



+ 2(D3) 



'j>cs p + 1 



pcs 



K 



m ° ^pcsfe. :I’) e igu 
y igu ° y 



+ S / %, \ © 

pcs(e ;rj rgu 

r rgu* & y 



2 (DR) a . 

+ yes p + 1 



i C 3 



yes 



I- = o 

li. 

z 



ycsie^Kfizu^ 



+ o / v0 e 

ycs(e ;i\; rgu 
rgu z 



Where : 



°(rcs) (e. ;!•:) d S (rcs)(e. ;T) 



igu' 



1£U' 



a (y»)(. Ulii K) - t 3 s (ycs)(e 

The remaining sensitivities have been s ini In rlv redefine.!. 

In writing equations to describe the gyro units, the gyre 
uncertainties are neglected due tc the relative!;, short system 
operatir g time. The rate gyr o units are assured to be r rorortior.a i 
mechanisms so that: 

(3-8) 

e - S X v; 

rgu rgu(/.je) x 

•A. 

' « S A, C ; 

rgu rgu(V.';e) y 



e = S 
r £U z 



rgu(V.';e ) W z 



where the sensitivity is the rate gyro sensitivity from input a*is 
angular velocity to output signal. 

The performance equations for the integrating gyro is developed 
in Derivation Summary 2 of Reference 2 . In terms of gyro input, 
output and spin reference a>es, the torque summation or the gyro 
girnbal is: 
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(B-V) 



• • 



”^( gim) A (gim) + A (I-B)CA “ A (vd)(A;K) A (gin.) 

• • 

+H A (I-B)(IA) C03 A (ginO " H A (l-B)(5aA) Rin A (gin) " ° 

The orientation of the gyro3 with resect to body axes is given 
in the following Table: 



Gyro 


1 Input Axis 


! Spin Reference Axis 


^ Output Pxi 


Roll 


X 

1 


i 

Y 

| 


z 


Pitch 


1 

V 


X 


i - z 


Yaw 


i z 


1 

V 


-X 



Table 5-1 



Gyro Orientations 



where the entries in the table refer to the missile body axes. 

Tne relationship between gyro gimhnl angle and gyro output 
signal is: 

(3-10) 

6 igu ~ S (igu)(A . ;e) A gira 

glill 

(3-9) and (3-10) xay be co.bined tc give the equation for gyro 
output signal. Scutions may be written for each axis by replacing 
the subscripts refe ring tc gyro ax°s by the appropriate body axis 
subscripts frcrr the Tr’ le ( 5-1) . Then, after rearranging (3-9) 
becorr.es : 

(3-11) 

*• • 

I (gim) e (igu) + * > (vd)( , .’. , ;K) e (igu) + K “y ®(igu) 



= H S/„ 



W - I 



(5g)(A;e) y “ (gira) (sg)(A;e) z 
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'(£im) e (igu) + J (vd)(v:.M)*(igu) + H W x e (igu) 



H °(Sg)(A;e) V 'y + I (gLn) S (Sg)(A;e) 



(gi.Ti) w (igu) (vd)(V.';K) (igu) 



’(igu). 



> (Sg)(A;e)"z A (gim) a (S S )(A;e) A x 
The equations (’-1), (3-2), (2-7) =n.1 (h-ll) cc:. 4 l^t-?] y c- a :ribe the 
ferfcrmance of the ; .royosed system. These can be sirylified by 
redefining certain of tee constants. 

Define : 

(H-l.) 

I - I 

(IE) = 

x I 

x 

I - I 

(iR) v = -T”" 

7 

I - I 

< IS >* = -V- 2 



/• _ s o let: 
(-5-13) 



‘(ir.t). 



= A/ . 



K 



m 



( in t } 

x ; 



= A 



m 



x 



K (in A ,) y 

I" “■ = A (int) 



‘(int) z 

T = A (int) 

z z 



7 



K 



n 



ir. 



= A 



7 



(int). 



The gyro angular moment irr* and the viscous darner sensitivity 
are equal for the gyro being used so that: 

(£-14) 



(CT), 



I . 
jSiE 



I . 



igU S (vd)(W;M) H 
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Substitution of (b-12) (B-13) (B-14) into the syst em equations has 
the effect of lunging the system sensitivity for each axis in the 
torque generating system equations, except for the signal generator 
sensitivity* If the gyro equations are expressed in terms cf gyro 
ginibal angle instead of voltage this sensitivity is also carried 
uu into the torque goner ting system The performance eqoati.ns 
for the systems then become: 



(5-1?) 



W - (IH) W • A/, .v + A 
x x y z (mt) 



v - (in) v: v. = a / . % + a 

y y x z lint} m 

y y 



Vi — ( J. . «. ) . . r . . - :i / . , \ + A 

z z x y mnt} m 

x 7 z z 



V. ^(rcs)(A . A ) H gim + ^rcsO'jA) A x 
z gun; ° x * 



2(M) 



DC S 



P r 1 



CCS 



OC3 



/\ = O 

m 



pcs (A . ;A ) M ri;r. 
r gun w y 



+ S 



> pcs( T, . r ; A ) *’y 



2(DR)^ 



4 . 1 



yes 



VC 5 



J ycs(A . ; A ) Vi.-r. 
“ v gun 7 -* ; 



+ i 



pcs(V.';A) "z 



(CT). 


A . + 


A 

n 


- a . 


- (ct). v; 


igu 


gim x 


gun 
^ x 


y 2 M X x 


1 ?U Z 


(CT). 


• • 

A . + 


A 

rt . 


+ v; a . = v.’ 


• 

+ (CT). 




gi~ 

y 


gl: V 


x gut y 


1C 11 2 


(CT). 


A . + 


A 

n 


v: a . = w 


+ (CT). 


lgVi 


gl>r. 2 




y KT-ta z 

J ^ z 


igu X 


: S 


(rc3)(A 


, ;a) ■ 




E . S ("il)(A;0 
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^(rcs)(v;;A) 




and identical definitions for the sensitivities held for the remaining 
axes of the system. 

The equations (B-15) are the set that wa solved on the digital 
computer in the present analysis. The values chosen for the various 
constants in these equ-tisns are given L: Table 3-1. 

Reduced Dynamics Case 

The simplified system which was iiscussei in Chatter ITT can be 
obtained from the equations (l-lr) b\ eliminating the dynamics and 
the non-linearities fro::, the last six equations. equators 



can then be cor.. 




( B-l;: ) 




X 






• • 
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The set of equations (B-16) was the set solved for the system with 



simplified dynamics. The constants here identical with those for 
the complete system which are given in Table 3-l« 

Single Axis P Linearized -Stability Analysis 

The question of system stability was approached by making a 
root locus plot of the linearized performance equations for the 
yaw axis. The yaw axis equations may be written: 

(B-17) 



i. - .925 v; v; + A / . s + A 

n x y int m 

• 'z z 



p 

1,25 



r 



p 

I 

I >>> 



4 1 



r 

* s ,, *: J k, 

z. ycs(n h) ( 

z gim 



ycs(V. ;A) 

~*^z S ycs(A . A) 

flT 

t 1 ■ 



D 






A . 

gxm 



= V ’z 4 33 ? 



..here numbers h ve been substitited from Table 3-1 f^r all exce t th 
sensitivities. These may be combined t:> give a single expression 
provided that the character of the r. on-linearities is preserved. 

That is, the order of fact rs in the non-linoar expression must b - * 
maintained. Then by substituting from the first ar..1 third equation 
into the second and their results: 

( B-1S) 



2 

P 


P P 




.. 1 

j 


_333 + p + \ 


n — + - +1 

|S25 25 




w - .925.*: v. - a, . . » 

z x y (mt) 1 



v; 



(pcs)(A gic ;A) j 333 + 



333 



1 * - WS5 (333 * p * '-yj 



v; 



where the ratio: 



(B-19) 

S 



cs(>.’;A) 
Vcs ( A ;A) 



- .1955 



is arbitrarily chosen 

If it is now assuned that there is nc n.oti n in roll and pitch 

ther^ 

(b-20) 

= 0 

x y 

T:ie equation (3-3?) then be c :,es linen r the performance function 
for tne yaw a/is *r.ay be v. ritten: 

O'— ai) 



IH-,,; 



jCs(A int ; ; 'x) 



2 

d 



333 



P / ">cr 4 



, 1 
* J j 



— 






. - o 


— 




2 

P 


333 


+ 1 


1 

\'f*‘ 


♦5s* i 


+ £ 


and trie chara 


cteristic 


equation 


is: 


(H-22) 












"1 


ycs( 




5 x . h->3 n 2 + . 




1 - 


2 r 






“ 2 






P 1 


.oo3p + l 


■J2- + 2_ 


4- 1 




L 




- 


^25 25 





y=s(A ; A ) 333 



top- -> 

• d*~ ,lQ?Fp * i 



= 1 - OF). 



= 0 



open loop 

The complex frequency plot of the open lco; performance fur.cticn 

(PF) q ^ as a function of S r cg /^ ,^\ is given in Fig. 3-3. The 
° /L3K gim’ ' 

plots and zeros of this plot are given below. 
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Foies 



Source 



P x - -333*3 
P 2’ P 3 = 0 



P,,F r = -12.5 1 j 21.65 



4’* 5 



Zeros 



Z 1 = - 32?. 9 



Z 2 = 5.1 



Integrating Gyro 

Integr ting Gyro an 1 Missile Equation 
Yav; Torcue Generator g System 



Rate Feedback 
Rate Fevdbrck 



As the open loop sensitivity is ir.cr*as'*J fron. zero it will 

be observed that the pole pair at the .rgin . . ve to the left while 

the pole r a i r du« to the torz^i gen -r.it ‘ t.; sterr nave to tv 

right. l:.e choice of tie o r er. lov s *; ni ti vitv S , > ;;*s 

ycsu . ;a; 

o ^ U1 

made so as to give both of th^se oscillatory modes the same damping 
r~tio, This n ears that both p'les a -at li^ . n the same radial 
line from the origin. By trial and* error t.nis radial was d^ter ined 
to oe the line for a dairying rati ^ of 0.14. Vfr.en both pole pairs 
lie or. this line the open lco; sensitivity v;as Jeter: *ro! as: 



(5-23) 




2 - . 0 



which is the value shown in Table 3-1* 
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